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ABSTRACT 

We present a first application of the recently proposed LITMUS test for magnetic helicity, as well as a thorough study of its applicability 
under different circumstances. In order to apply this test to the galactic magnetic field, the newly developed critical filter formalism 
is used to produce an all-sky map of the Faraday depth. The test does not detect helicity in the galactic magnetic field. To understand 
the significance of this finding, we made an applicability study, showing that a definite conclusion about the absence of magnetic 
helicity in the galactic field has not yet been reached. This study is conducted by applying the test to simulated observational data. 
We consider simulations in a flat sky approximation and all-sky simulations, both with assumptions of constant electron densities and 
realistic distributions of thermal and cosmic ray electrons. Our results suggest that the LITMUS test does indeed perform very well in 
cases where constant electron densities can be assumed, both in the flat-sky limit and in the galactic setting. Non-trivial distributions 
of thermal and cosmic ray electrons, however, may complicate the scenario to the point where helicity in the magnetic field can escape 
detection. 
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1. Introduction 



ducting plasma becomes 



Helicity is of utmost interest in the study of astrophysical mag- 
netism. Mean field theories for turbulent dynamos operating 
in the galactic interstellar medium have been successful in ex- 
plaining how the observed magnetic field strengths are main- 
tained (e.g. Subramanian 2002). These theories predict that he- 
licity is present on small scales in interstellar magnetic fields. 
Observationally detecting or excluding helicity in these fields 
would therefore either strongly suggest that these theories are 
valid or indicate that there are some flaws in them. 

However, since helicity is a quantity that describes the three- 
dimensional structure of a magnetic field and most observation 
techniques produce at best two-dimensional images leading to 
an informational deficit, it has thus far largely eluded observers. 
Previous work on the detection of magnetic helicity in astro- 
physical contexts has focused mainly on either magnetic fields 
of specific objects, such as the Sun (see e.g. |Zhang||2010| and 
references therein) or astrophysical jets (cf. e.g. EnBlin 2003 ; 
Gabuzda et al. 2004), or cosmological primordial magnetic fields 



B j oc B (V x B) . 



(1) 



(eTgTTKahniash vili & Ratra|2005||K ahniash vili et al.|2 005). Two 



exceptions are the work by Volegova & Stepanov (2010)1, in 
which the use of Faraday rotation and synchrotron radiation for 
detecting magnetic helicity was suggested for the first time, and 
the work of Kahniashvili & Vachaspati (2006), in which the use 
of charged ultra high energy cosmic rays of known sources is 
suggested for probing the three-dimensional structure of mag- 
netic fields through which they pass. However, the sources of 
ultra high energy cosmic rays are not known yet and the appli- 
cability of this test is therefore limited. 

The LITMUS (Local Inference Test for Magnetic fields 
which Uncovers heliceS) procedure for the detection of mag- 
netic helicity suggested by Junklewitz & EnBlin ( 2010[ ) probes 
the local current helicity density B ■ j, which for an ideally con- 



Here, the magnetic field is denoted by B and the electric cur- 
rent density by /. The test uses measurements of the Faraday 
depth and of the polarization direction of synchrotron radiation 
to probe the magnetic field components along the line of sight 
and perpendicular to it, respectively. Its simple geometrical mo- 
tivation should make it applicable in a general setting, provided 
these quantities can be measured. The results depend only on 
the properties of the magnetic field along a line of sight and are 
therefore purely local in the two-dimensional sky projection. Our 
aim is to test this idea on observational as well as on simulated 
data, thereby determining the conditions under which the test 
will yield useful results. 

This paper is organized as follows. In Sect.|2j the basic equa- 
tions used in the LITMUS test are reviewed. They are applied 
to observational data describing the galactic magnetic field in 
Sect. [3] with special emphasis on a sophisticated reconstruction 
of the Faraday depth, described in Sect. 3.2 Section |4] is de- 
voted to a thorough general assessment of the test's reliability. 
To this end it is applied to simulated observations of increasing 



complexity. Section 4.1 describes the application in a fiat sky 
approximation, whereas Sect. 4.2 examines all-sky simulations, 
finally arriving at complete simulations of the galactic setting in 
Sect. |4.2.2[ where realistic electron distributions are added. We 
discuss our results and conclude in Sect. [5] 



2. The helicity test 

For a thorough introduction into the ideas behind the LITMUS 
test, the reader is referred to Jun klewitz & E nBlin (2010j). Here, 
we only summarize the resulting equations. 
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On the one hand side, synchrotron emission produced by 
cosmic ray electrons is used to probe the magnetic field com- 
ponent perpendicular to the line of sight. Its polarization is de- 
scribed by the complex field 



3.1. Observational data 



p = Q + iU = \P\e 



2<> 



(2) 



where Q and U are the usual Stokes parameters quantifying the 
linearly polarized components of the radiation with respect to 
some orthogonal coordinate system and^ is the polarization an- 
gle with respect to the first coordinate direction. On the other 
hand, the Faraday depth 



0, 



J LOS 



n e B ■ dl 



(3) 



is used to probe the magnetic field component parallel to the line 
of sight (LOS). 

A helical magnetic field will lead to a gradient of the Faraday 
depth that is parallel to the polarization direction of the syn- 
chrotron emission, as was argued in Junklewitz & EnBlin (2010). 
In order to compare the directions of the two quantities, this gra- 
dient is also formulated as a complex field 



G = ((Vtp) x + i(Vtp) y f = \G\, 



,2/ft' 



with 



a = arctan 



cm 



(4) 



(5) 



where the indices x and y denote its components with respect to 
the coordinates used. The helicity test that is performed in this 
work consists simply of multiplying G with the complex conju- 
gate of the polarization P*. If the two angles x an d ot differ by a 
multiple of n (i.e. the gradient and the polarization direction are 
parallel), the product will be real and positive. If they differ by 
an odd multiple of n/2 (i.e. the two directions are perpendicu- 
lar), it will be real and negative. Any orientation in between will 
produce varying real and imaginary parts in the product. Thus, 
observational directions along which a magnetic field is helical 
are indicated by a positive real part and a vanishing imaginary 
part of the product. Averaging over all observational directions 
will give an indication of the global helicity of the field. 

It was furthermore shown by Junklewitz & EnBlin ( 2010| > 
that the ensemble average of this product over all magnetic field 
realizations given a magnetic correlation tensor (and therefore a 
helicity power spectrum) is a measure for the squared integrated 
spectral current helicity density 



<GP*> £ 



'd*«f) 



(6) 



with large scales weighted more strongly than small scales. 



3. Application to galactic observations 

In this section, we try to answer the question whether the mag- 
netic field of the Milky Way is helical by applying the LITMUS 
test to the available observational data. Since the magnetic field 
is localized in a region that surrounds the observer, all relevant 
quantities will be given as fields on the sphere S 2 , i.e. as func- 
tions of the observational direction, specified by two angles i? 
and <p, which are taken to represent the standard spherical polar 
coordinates in a galactic coordinate system. 



For the synchrotron emission, we use the data gathered by the 
WMAP satellite after seven years of observation^] described in 



Page et al. 



( 2007| >. Since the foreground synchrotron emission 
is most intense at low frequencies, we use the measurement in 
the K-Band, which is centered at a frequency of v = 23 GHz. 
Furthermore, we assume that the detected polarized intensity is 
solely due to galactic synchrotron emission. Thus, the Stokes 
Q and U parameter maps (defined with respect to the spherical 
polar coordinate directions <?,<> and in the galactic coordinate 
system) can be simply combined according to Eq. |2]l to give the 
complex quantity P whose argument is twice the rotation angle 
of the plane of polarization with respect to the e^-direction 



X{&, <f) = - arctan 



/Im(P(#, ip)) 
\Re(P(#, ip)) 



(7) 



(cf. |Junklewitz & EnBlin(20T0) l. 

There are several depolarizing effects that have to be con- 
sidered when dealing with polarization data. Faraday depo- 
larization, which is important only at low frequencies due to 
the proportionality of the Faraday rotation angle to the square 
of the wavelength, can be safely neglected in the K-Band. 
Depolarization effects due to different magnetic field orienta- 
tions along the line of sight are certainly present. However, 
they are present as well in the numerical test cases presented 
in Sect.|4] which yield good results. Additionally, depolarization 
due to the finite beam-size of the WMAP satellite and the finite 
pixel size of the polarization maps used in this study can play 
a role. The only way to limit this effect is to use higher resolu- 
tion maps, ultimately necessitating the use of Planck data in the 
future. 

In order to construct a map of the Faraday depth, we use 
the cat alog of rotation measurements provided by |Taylor et al.| 
(20090 These provide an observational estimate of the Faraday 
depth for certain directions in the sky where polarized radio 
point-sources could be observed. Since the catalog encompasses 
a large number (37 543) of point-sources, it paints a rather clear 
picture of the structure of the Faraday depth. However, Earth's 
shadow prevents observations in a considerably large region 
within the southern hemisphere. 

3.2. Reconstructing the Faraday depth map 

The reconstruction is conducted according to the critical filter 
method first presented in|Enfilin & Frommert (20101. A more 



elegant derivation of the same filter can be found in EnBlin & 
Weig ( 2010] >. Since this formalism takes into account available 
information on the statistical properties of the signal in the form 
of the power spectrum, it is able to interpolate into regions where 
no direct information on the signal is provided by the data, such 
as the shadow of Earth in this case. Furthermore, it takes into 
account the available information on the uncertainty of the mea- 
surements. All in all it is expected to lead to a reconstructed map 
of the Faraday depth that is much closer to reality than e.g. a 
map in which the data were simply smoothed to cover the sphere. 
Small-scale features that are lost in such a smoothing process are 
for example reproduced by the critical filter algorithm. 



1 The data are available from NASA's Legacy Archive for Microwave 
Background Data Analysis at http : //lambda . gsf c . nasa . gov. 

2 The catalog is available at 
http : //www . ucalgary . ca/ras/ rmcatatlogue. 
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Fig. 1. Vertical galactic profile p(&) of the Faraday depth. 



3.2.1. Data model 

The field that is to be reconstructed here is the sky-map of the 
Faraday depth. In order to apply the critical filter formula, the 
signal should be an isotropic Gaussian field. Since the Faraday 
depth clearly is larger along directions passing through the galac- 
tic plane, the condition of isotropy is not satisfied. Therefore, a 
vertical profile is calculated by binning the observations into in- 
tervals [§i, §i + A§), calculating the root mean square rotation 
measure value for each bin and smoothing the resulting values 
to obtain a smooth function p(§). The result is shown in Fig.[T] 
This profile is used to approximatively correct the anisotropies 
induced by the galactic structure and the resulting signal field 



(fi) - 



P {§) 



(8) 



is assumed to be isotropic and Gaussian with a covariance ma- 
trix S . The Gaussian covariance matrix is determined solely by 
the angular power spectrum coefficients Q, the reconstruction of 
which is part of the problem at hand. 

The data d, i.e. the rotation measure values in the catalog, 
are taken to arise from the signal s by multiplication with a re- 
sponse matrix R, which consists of a part encoding the specific 
directions in which the signal field is probed in order to pro- 
duce the measurements and another part that is a simple multi- 
plication with the vertical profile /?(#). Additionally, a Gaussian 
noise component n is assumed with a covariance matrix N = 
diag(<x 2 , cr\, . . . ), where cr, is the one sigma error bar for the ith 
measurement in the catalog. Thus, the data are given 

d = Rs + n - Rps + n. (9) 

Recent discussions in the literature (see e.g. Stil et al. 201 1 1 have 
shown, however, that the error estimates as quoted in the cat- 
alog are probably too low. In addition, any contribution to the 
measured data from intrinsic Faraday rotation within the sources 
will further increase the error budget since the signal field in this 
context is only the contribution of the Milky Way. We therefore 
adapt the error bars of Taylor et al. (2009) according to the for- 
mula 



^(corrected) = ^f^f + ^mt)f 



(10) 



3 The discretized version used in the implementations is d, = 
YjjRijSj + n i = YijRijPj s j + n h where the index j determines a pixel 
on the sphere, so that Sj = s(&j, <pj) and pj = p(&j). 



respectively. Here, the factor fa- accounts for the general under- 
estimation of the errors in the catalog of |Taylor et al. (2009), 
whereas the additive constant cr (lnt) represents the average con- 
tribution of the sources' intrinsic Faraday rotation. As numeri- 
cal values, we use f a = 1.22, which was found by Stil et al. 
((20TT) by comparing the data of Taylor et al. ( 2009 1 and Mao 



et al.| ( 20T0] l on sources that are contained in both catalogs, and 
cr 1 " 11 - 1 = 6.6 m" 2 , which corresponds to the upper end of the num- 
bers found by Schnitzeler ( 2010| l. Since the error contribution 
from the internal Faraday rotation is not correlated with the mea- 
surement error, the two contributions add up quadratically in the 
noise covariance matrix. 



3.2.2. Reconstruction method 

For signal and noise fields with a Gaussian prior probability dis- 
tribution, the posterior probability distribution of the signal field 
is again a Gaussian, i.e. it is of the form !P( sic/, C/) = Q(s-m,D), 
where a multivariate Gaussian probability distribution with co- 
variance matrix X is denoted by 



1 



\2nX\ 



1/2 



exp 



1 



(11) 



where the t symbol denotes a transposed and complex con- 
jugated quantity. In order to reconstruct the mean signal field 
m = (s), where the brackets denote the posterior mean, it is also 
necessary to reconstruct the angular power spectrum C\. To do 
so, the critical filter formulas 



m — Dj 



(12) 



and 



C/= 27TT tr (( mm1 +D ) S >) (13) 

are iterated, starting with some initial guess for the power spec- 
trum. Here, the signal covariance matrix is expanded as S = 
2/ C[S i, where Si is the projection onto the spherical harmonic 
components with index /. Furthermore, D is the posterior covari- 
ance matrix, 

D = (S~ l +R t N~ l Ry 1 , (14) 
and j is the information source term, 



R^N^d. 



(15) 



Since the critical filter is on the brink of exhibiting a per- 
ception threshold (cf. EnBlin & Frommert 2010) and it is gen- 
erally more desirable to overestimate a power spectrum enter- 
ing a filter than to underestimate it, the coefficients Q are sub- 
jected to a procedure in which the value of C; is replaced by 
max{C/_i,C/, C/ + i) after each iteration step. The advantage of 
overestimating the power spectrum can be seen by considering 
the limit of C/ — > oo in Eq. (12\ and ( |14) , For high values of C/, 
the first term in Eq. ( fl4| > can be neglected and Eq. ( fT2) becomes 



mci- 



(16) 



Thus, by overestimating the power spectrum the importance of 
its exact shape is diminished and the reconstruction will instead 
follow the information given directly by the data more closely. 
Considering an extreme underestimation of the power spectrum, 
C/ — > 0, on the other hand, would lead to 



m c ,^Q = 0, 

suppressing the information given by the data. 



(17) 
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Fig. 2. Results of the reconstruction of the signal field and Faraday depth. The left column shows the posterior mean of the signal 
field m (panel (a)) and its one-sigma uncertainty (panel (c)). The right column shows the resulting map of the Faraday depth 

pm (panel(b)) and the corresponding one-sigma uncertainty Jp 2 D (panel (d)) in irT 2 . 




0.0001 - 



Fig. 3. Angular power spectrum of the signal field s. The solid 
curve shows the resulting power spectrum as calculated with 
the critical filter formalism. The dashed line depicts the signal 
power spectrum used to generate large-scale uncertainty correc- 



tions for the calculation of G according to Sect. 3.2.3 



3.2.3. Calculating {G) g(s _ 



m,D) 



Once the posterior mean of the signal is reconstructed, the corre- 
sponding field G can be calculated according to Eq. (|4]), assum- 
ing (p - m. However, since Eq. (|4]> is a nonlinear transformation 
of the Faraday depth, that is not the same as the posterior mean 



of G. This can be estimated by averaging over a large number of 
samples, according to 



(G(s)) gt 



1 



(s-m,D) 



samples 'j^y 



>l11I1] Ml. > 



(18) 



Here, the samples s% are drawn from the posterior probabil- 
ity distribution Q(s — m,D). In our implementation we use the 
method of Jasc he etaL| (|2010) to generate these samples, i.e. we 
generate corrections y k to the posterior mean m by drawing sig- 
nal and noise realizations from the respective prior probability 
distributions, 

h*^Pis\S), (19) 
n k ^ P(n\N), (20) 
combining them into a data realization according to 

d k =Rs k + n k , (21) 

and calculating the difference between the signal realization and 
the Wiener filtered data realization 



y k = s k - DR*N~ l d k . 
The signal realizations 

Sk = m + y k 



(22) 



(23) 



are shown by Jasche et al. (2010 1 to follow the posterior proba- 
bility distribution. For reasons of numerical feasibility, we con- 
sider only large scale corrections to the posterior mean in our 
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0.01 




Fig. 4. Results of the calculation of the gradient field G and the application of the LITMUS test. The first row shows estimates for 
the ablolute values of the field G. Panel (a) shows G as calculated directly from the posterior mean m of the signal field, whereas 
panel (b) takes into account large-scale uncertainty corrections as described in Sect. 3.2.3 Panel (c) shows the absolute value of the 



polarized intensity P as measured by the WMAP satellite and panel (d) shows the result of the LITMUS test, namely Re {GP*), in 
arbitrary units, using G and P as shown in panels (b) and (c), respectively. Note the logarithmic color code of panels (a)-(c). 



calculation, i.e. we introduce a small-scale power cut-off for the 
signal power spectrum that describes the prior probability distri- 
bution used in Eq. ( |19) , 

Thus, given m and D, the posterior mean for the Faraday 

depth is given by (cf>) = pm, its one sigma error bars by ±p V5, 
and the posterior mean for the field G by the above sampling ap- 
proximation, where the signal samples are multiplied with the 
galactic profile function p{9) to give realizations of the Farady 
depth. Here, D is the vector that contains the diagonal elements 
of the matrix D. 

3.3. Results 

Figures [2] - 14] summarize the results of the Faraday depth re- 
construction and the application of the LITMUS test to these 
data. All calculations are conducted at a HEALPdQ resolution 
Nside = 64. The left column of Fig. [2] shows the reconstructed 
signal field as well as its uncertainty as they are calculated using 
the critical filter method. Evidently, the reconstruction method 
is able to extrapolate from the available information into regions 
where no data are taken, i.e. the Earth's shadow in the lower right 
of the projection. However, only structures on scales comparable 
to the extent of the region without information are reconstructed 
within it and the reconstruction's uncertainty becomes large in 
this region, as well as near the galactic poles, where the signal- 
response to noise ratio of the data is low. 

4 The HEALPix package is available from 
http : / /healpix . j pi . nasa . gov. 



The outcome for the angular power spectrum is shown in 
Fig. [3] Clearly, the critical filter predicts a large amount of power 
on small scales. Therefore the rather high values seen in Fig.|2]c 
are mainly due to uncertainty on small scales. This may in part 
be due to an underestimation of the error bars of extreme data 
points which are not accounted for by our rather crude correction 
of the error bars as described by Eq. (fT0])p| 

The right column of Fig.[2]shows the reconstructed Faraday 
depth and its uncertainty, obtained from the left column of Fig. [2] 
by a simple multiplication with the galactic profile function 
p{&). Clearly, the field {(/>) takes on only low values within the 
information-less region, whereas, again, its uncertainty is espe- 
cially large there. 

Figure |4] shows the two fields used in the LITMUS test as 
well as its result. The field G calculated according to the sam- 
pling prescription, Eq. ( fT8] >, is shown in panel (b). In the sam- 
pling procedure, large-scale uncertainty corrections that are de- 
scribed by the signal power spectrum shown in Fig. [3] are taken 
into account. For comparison, the field G as calculated from the 
posterior mean of the signal field without any uncertainty cor- 
rections is shown in panel (a). The main effect of the uncertainty 
corrections is to increase the values within the region where no 
data points were measured and in the polar regions. The spatial 
structure of G, however, remains largely unchanged. 



5 A thorough study of the problem of reconstructing a signal field 
with unknown power spectrum from data with unknown error bars will 
be presented in a future paper. 
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field without divergence, its frequency components parallel to 
the respective ^-vector are then subtracted 



Fig. 5. Rotational curve for Re (GP*)^ in arbitrary units. 



Finally, the real part of the product GP* is shown in the 
last panel of Fig. [4] By visual inspection, it is not immediately 
clear whether positive values dominate this image. Taking the 
spatial average Re(GP*) s i does in fact yield a negative value. 
Nevertheless, we performed another test for helicity. In order to 
assess helicity on small scales, the two fields G and P are rotated 
with respect to one another about an angle /3 around the galactic 
axis and their product is again averaged over the whole sky. For 
a helical magnetic field, this procedure should result in a curve 
Re (GP*) s i (J3) that exhibits a maximum at [3 — 0, since the cor- 
relations between the fields G and P are expected to be local. The 
resulting curve for this test is shown in Fig. [5] Not only does it 
not take on its global maximum at /3 — 0, but the point /3 = 
does not seem to be special in any way. Therefore, no indication 
of helicity is found. 

We also performed the analysis with a different function 
p("ff), created by doubling the smoothing length in its calcula- 
tion. The results of the reconstructed map of the Faraday depth 
changed only slightly, i.e. the variance of the difference between 
the reconstructions with the two different profiles is 0.65% of the 
variance of the original (^)-map shown in Fig.[2]b. 

4. Application to simulated data 

In order to check whether the nondetection of helicity in the pre- 
vious section allows the conclusion that the galactic magnetic 
field is in fact nonhelical, we now apply the same helicity test 
to a number of artificially generated magnetic fields with known 
helicity. 

4.1. Planar implementation 

The most simple setting that can be considered is the observation 
of a well localized magnetic field structure. In the limit of great 
distances between the magnetic field under consideration and the 
observer, the lines of sight penetrating the field become parallel. 
We assume the field to be contained in a cubic box which is 
oriented along the lines of sight. 

The field in the box is generated by the garfields code (first 
applied in Kitaura & EnBlin 2008j). This code draws the three 
cartesian components of the magnetic field in Fourier space in- 
dependently of a common power spectrum, assumed here to be 



#div-free(*) K B ( k ) ~ k 



k ■ B{k) 

k 2 ' 



(24) 



A degree of helicity is then imprinted onto the field by applying 
the formula 

d ti\ d /i\ , X ^div-free(*) 

»div-free,helW K "div-freeW + V ; . 

K 

where r\ — leaves the field unaffected and rj — ± 1 produces the 
highest degree of helicity. 

Finally, we assume the thermal and cosmic ray electron den- 
sities to be constant throughout the box. Thus, the observables 
Q, U, and (f> can be obtained by simply integrating the appro- 
priate magnetic field components along the box direction asso- 
ciated with the line of sight. Then the complex quantities G and 
P are easily calculated and multiplied, yielding two-dimensional 
images of GP* . This procedure is conducted for various realiza- 
tions of random magnetic fields both without helicity (77 = 0) 
and with maximal helicity (77 = 1). 



4.1.1. Results 

The resulting images for one random magnetic field realization 
are shown in Fig. [6] The cube was discretized for the calcula- 
tion into 512 3 pixels. It can already be seen by eye that positive 
values of Re(GP*) dominate in the case with maximal helic- 
ity (panel (a)), whereas in the case without helicity (panel (b)), 
positive and negative values seem to be roughly equally repre- 
sented. Calculating the spatial averages over the whole square 
yields Re<GF% = 2.3 • 1(T 3 and Re<GP*) D = -1.5 ■ 1(T 4 for 
the case with and without helicity, respectively. 

We calculated this spatial average for the results of the 
LITMUS test applied to 100 different random magnetic field re- 
alizations, both with and without helicity. Averaging these val- 
ues for the helical fields and for the nonhelical fields separately 
yields a positive value in the helical case. Normalizing all values 
such that this average becomes equal to one yields 

<Re<GP*) D ) samples = 1.0, o- Re(cnn =0.63 

in the case with helicity and 

<Re<GP*) D ) samples = -0.025, cr Re(cnn = 0.34 

in the case without. Clearly, the LITMUS test yields positive 
results if applied to helical fields, whereas its results fluctuate 
around zero if applied to nonhelical fields. This is exactly the 
behavior that should be expected and the basic functioning of 
the LITMUS test is thereby demonstrated in this setting. 

4.2. Spherical implementation 

As a next step, the applicability of the LITMUS test is checked 
for magnetic fields surrounding the observer, as in the case of 
the galactic field. Again, several sets of mock observations are 
produced. These simulations are conducted using the hammurabi 
code (see Wae lkens et al.|2009] l in connection with the garfields 
coder] The hammurabi code allows for a large scale analytic 
field model and an additional Gaussian random field component, 



-5/3-2 



a Kolmogorov-type spectrum of the form Ps{k) k 
cording to Gaussian statistics. In order to produce a magnetic 



Both codes are available from 

http : //www.mpa-garching .mpg . de/hammurabi/hammurabill. 




Fig. 6. Maps of Re (GP*) for a particular magnetic field realization in arbitrary units. Panel (a) shows the case with maximal helicity 
(77 = 1), panel (b) the one without helicity (77 = 0). 



which can be generated by the garfields code with a preset de- 
gree of helicity as described in Sect. |4. 1| It then integrates the 
different field components, weighted with the appropriate elec- 
tron density, along radial lines of sight and produces sky maps 
of simulated observations of the Stokes parameters Q and U and 
Faraday depth <p (among others), thus providing all necessary 
ingredients to perform the LITMUS procedure. 



4.2.1. Constant electron densities 

For simplicity, we start again by setting the densities of the ther- 
mal electrons and the cosmic ray electrons to constant values 
throughout the simulated galaxy. 



Gaussian random field. First, in order to apply the test to a 
field with a well-defined degree of helicity, the field strength of 
the large scale analytic component is set to zero, such that the 
simulated galactic field is a purely random one with a chosen 
degree of helicity. As in the planar case, we choose either no 
helicity (77 = 0) or maximal helicity (rj = 1). 

Figure [7] shows the maps of Re (GP*) for one particular 
Gaussian random magnetic field with a power law index of -5/3, 
with helicity and without helicity, respectively. As for the planar 
implementation, one can immediately see that positive values 
dominate in the helical case and positive and negative values are 
roughly equally represented in the nonhelical case. 

Figure [8j shows the results of the rotational test described in 
Sect. |3.3| for this particular magnetic field realization. Clearly, 
the spatial average Re(Gf*)^2 takes on a sharp maximum at 
/3 = and is positive in the case with helicity, while it does 
not have a maximum there and in fact happens to be negative in 
the case without helicity. This is exactly the result expected from 
the LITMUS test. 



The sharpness of the peak in Fig. [Syndicates that the helicity 
is to be found in small-scale features. As a test of this assertion, 
we calculated the spherical multipole components G/„, and P/„, 
of the complex gradient and polarization fields. Note that the 
spatial average over the product of the fields is proportional to 
the sum of the products of the multipole components, i.e. 



Re(GP*)« = —Re 
4tt 



\ 1=0 m=-l 



(26) 



where / max is determined by the finite resolution of the map. If 
we now neglect the first terms in the sum, i.e. the small-Z contri- 
butions, we arrive at a spatial average over the product in which 
all large-scale features were neglected. The resulting rotational 
curves for the same magnetic fields used for Fig.[8jare shown in 
Fig. [9] where only multipole moments with / > 25 were consid- 
ered. It can be seen that this procedure further sharpens the peak 
at /3 = and strengthens it relative to other local maxima in the 
curves, thus facilitating the detection of small-scale helicity. The 
same result for the observational data studied in Sect.[3]is shown 



in Fig. 10 Clearly, there is still no sign of helicity in this case. 

Furthermore, we created a set of 100 different Gaussian ran- 
dom magnetic fields, performed the LITMUS test, and calculated 
the spatial average Re (GP*) $i for all of them. Each field realiza- 
tion was considered in a version without helicity and a version 
with maximal helicity, just as in the case of the planar imple- 
mentation. Averaging over the 100 samples yields again a pos- 
itive value in the helical case. Normalizing all values such that 
this average is equal to one yields 

<Re(GP*> 52 ) samples = 1.0, cr Re(GP%2 = 0.25 

in the helical case and 

<Re<GP*>^) samples = -0.27, cr Re(GP%1 = 0.23 
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Fig. 7. Maps of Re (GP*) for a particular magnetic field realization in arbitrary units in a spherical setting. Panel (a) shows the case 
with maximal helicity (77 = 1), panel (b) the one without helicity (77 = 0). 




a. 




Fig. 8. Rotational curve for Re(Gf*)^2 in arbitrary units. The Fig. 9. Same as Fig. [i] only with large-scale contributions ne- 
solid and dashed lines depict the results for the same magnetic glected. 
field with helicity parameter r\ — 1 and r\ — respectively. 
Constant electron densities are assumed. 



in the nonhelical case. This clearly underlines the success of the 
LITMUS test in the spherical setting. 

Since the magnetic field of the Milky Way is known to be 
much stronger in the vicinity of the galactic plane than in the 
halo region, we repeated the analysis with random realizations of 
Gaussian magnetic fields that are confined to the region within a 
vertical distance of 8 kpc of the galactic plane. The results of the 
LITMUS procedure for these fields still turn out to be a reliable 
indicator for magnetic helicity. 



Large scale field models. As a next step, the (helical or non- 
helical) random magnetic field component is switched off com- 
pletely and replaced by an analytic large scale magnetic field 
model. Several sets of simulations are performed using differ- 
ent models for the galactic large scale field. For these analytic 
models, the current helicity can be calculated directly, giving an 
expectation as to whether the helicity test should produce posi- 
tive results or not. 



a, 

-wi 




Fig. 10. Same as Fig. [5] only with large-scale contributions ne- 
glected. 



( 12007) , 



i.e. 



The results of the rotational helicity test are shown in Fig. 1 1 
for the large scale magnetic field model described in Page et al. 



B(r, ip, z)=B [ cos (tff(r)) cos (x(z)) e r 
+ sin (i/f(r)) cos (x(z)) e v 
+ sin CvrCz)) e z ] 



(27) 
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S 0.6 P 




Fig. 11. Rotational curve for Re(Gf*)^2 in arbitrary units. The 
solid line corresponds to the best-fit magnetic field model of 
Page et al. ( 2007 1, the dashed line to a model where the parame- 
ter xo is set to zero. Constant electron densities are assumed. 



in galactic cylindrical coordinates, where 



and 



ij/(f) = i// + ifri In 



X(z) = Xo tanh 



8kpc 



lkpc 



(28) 



(29) 



The solid curve corresponds to the parameters favored by Page 
eFaLl ( |2007l >, namely xo = 25°, i/f = 27°, and if/ x = 0.9°. This 
corresponds to a simple fiat spiral in the galactic plane which 
becomes more and more screw-like with vertical distance z from 
the galactic plane, so that a slight degree of helicity is inherent 
in the field geometry. This can be verified by direct calculation 
according to Eq. ([TJ, yielding 



JB = Bl 



sinCr(z)) cos (x(z)) 



(sin(^(r)) + ^icos(^(r))), (30) 



which is nonzero for any generic point away from the galac- 
tic plane. The resulting line in Fig. 11 is not a clear indication 
for this helicity. However, the curve is nevertheless sensitive to 
the angle xo, which produces the helicity. Lowering its value, 
i.e. making the spirals more and more parallel to the galactic 
plane, changes the results of the LITMUS test. The extreme case 
of^o = 0, i.e. B z = everywhere, for which the value of Eq. ( f30| > 
becomes zero everywhere, is also shown in Fig. TT| This curve's 



value at B — is even more distinct from its maximum than in 
the case of the solid line. This example shows that while the re- 
sults react in a systematic way on changes in the parameters, the 
test is not suited to detect helicity on the largest scales. 

We used the model of |Page e t al. (2007) in the demonstra- 
tion of this effect mainly because of its mathematical simplicity. 
More sophisticated models can be found e.g. in Sun et al. (2008 ), 
|Jansson et al.| ( |200"9"| i, |Jaffe et al.| ( |2010| l, and references therein. 



4.2.2. The role of the electron densities 

In order to get closer to a realistic model of the Milky Way, as a 
next step we replace the constant electron densities with realis- 
tic models. The hammurabi code allows the use of the NE2001 
model for the thermal electron density (cf. |Cordes & Lazio 2002 



a. 



Fig. 12. Same as Fig. [i] only with realistic electron densities 
used in the calculation. 



6. 




Fig. 13. Same as Fig. [IT] only with realistic electron densities 
used in the calculation. 



|2003 1 to compute the Faraday depth and several analytic models 
for the cosmic ray electron density to compute the synchrotron 
emissivity (see also Waelkens et al. 2009). In the calculations 
performed to obtain the results presented here, the cosmic ray 
electron distribution model of Pag e et aL] ( |2007[ l was used. 

Gaussian random field. The resulting rotational curves, as cal- 
culated with the realistic electron distributions, for the case of 
the Gaussian random field are shown in Fig. 12 Neither in the 
case with helicity (77 = 1), nor in the case without helicity (77 = 0) 
does the curve take on its maximum at B — 0. This is the expected 
result in the latter case but contradicts the expectation in the for- 
mer one. Therefore, the helicity imprinted onto the small scale 
magnetic field clearly fails to be detected by the test applied. 

Large scale field models. The result for the large scale mag- 
netic fields is shown in Fig. 13 There actually seems to be a 

now. This is true, however, for 
0) as well as for the 



maximum in the vicinity of B 
the case of the planar spiral model (xo 
model with helicity (xq = 25°). Therefore, the proposed helicity 
test might under certain circumstances even indicate helicity on 
large scales where there is none, if the observer is surrounded by 
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the field. Other magnetic field models with planar spirals, such 
as the bisymmetric (i.e. B(r, <p, z) = -B(r, if + n, zj) spiral model 
of |Stanev| ( |1997| l lead to similar results. 



ing his visit to Germany, which was supported by the DFG-RFBR grant 08-02- 
92881. We thank the anonymous referee for constructive comments. 



5. Discussion and conclusion 

The present work presents the first application of the LITMUS 
test for magnetic helicity proposed by Junklewitz & EnBlin 
( |2010| l to actual data. The application of the test involves thor- 
oughly reconstructing a map of the Faraday depth distribution, 
calculating its transformed gradient field G, creating a map of 
GP*, averaging over this map, shifting the two fields with respect 
to each other to see whether any signal vanishes, and filtering 
out large-scale contributions for a better detection of small-scale 
helicity. This procedure, applied to observations of the Faraday 
depth and polarization properties of the synchrotron radiation 
within our own galaxy in Sect. [3] does not show any signs of 
helicity in the Milky Way's magnetic field. 

In order to assess the significance of this, the applicability of 
the test was probed in different artificial settings. The complex- 
ity of these settings was increased bit by bit to find out under 
what circumstances exactly the LITMUS test yields reliable re- 
sults. It was found that meaningful results can be achieved if 
the electron densities do not vary on the scales of the magnetic 
field, both in the regime of magnetic field structures whose dis- 
tance from the observer is much greater than their extension, as 
shown in Sect 



4.1 



and in the regime of magnetic fields sur- 
rounding the observer, as shown in Sect. |4.2.1| We showed that 
the performance of the LITMUS test with regard to small-scale 
helicity is further improved by dropping the first few terms in 
Eq. ( p6| ). However, indications of helicity on large scales are un- 
reliable, as shown in Sect. |4.2.1| Furthermore, it was demon- 
strated in Sect. 4.2.2 that any non-trivial electron density may 
distort the outcome of the test to a point where even small-scale 
helical structures fail to be detected. This is not too surprising 
since e.g. a variation in the thermal electron density will intro- 
duce a gradient in the Faraday depth that is not caused by the 
magnetic field structure. Therefore the nondetection of helicity 
for the galactic magnetic field does not necessarily mean that the 
field is nonhelical on small scales. It may be the case that small- 
scale fluctuations of the electron density introduce effects in the 
observational data that prevent the detection of helicity. 

So, as a natural next step, the hunt for helicity in astrophysi- 
cal magnetic fields should focus on a region that is small and/or 
homogeneus enough for the assumption of constant electron 
densities to hold at least approximatively. Although our work 
has shown that the helicity test that we studied is not suitable 
for all astrophysical settings, we are confident that it may nev- 
ertheless yield useful results if applied in a setting with constant 
electron densities. 

As a side effect of this paper, it was demonstrated in Sect. 3.2 
that the method proposed by EnB lin & Frommert] ( |2010| l to re- 
construct a Gaussian signal with unknown power spectrum is 
very well suited for practical application. 
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